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Abstract. 

We present an updated timing solution and an analysis of the profile evolution - including precession and beam shape 
- of the young, relativistic binary pulsar J1906+0746. The 144-ms pulsar, in a 3.98-hour orbit with eccentricity 0.085 [l|], 
was initially discovered during the early stages of the ALFA (Arecibo L-band Feed Array) pulsar survey |2] using the 305- 
metre Arecibo telescope and was subsequently found in archival Parkes Multibeam Survey data. We have since been regularly 
monitoring the system using the Arecibo and Green Bank telescopes, and include data from the Jodrell Bank, Parkes, Nancay 
and Westerbork telescopes. The nature of the binary companion will also be discussed based on improved estimates of the 
total and companion masses obtained from the updated timing solution. 
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GAUSSIAN FITTING AND PROFILE 
EVOLUTION 

In order to better understand the profile evolution and in 
turn remove some timing noise from the system, we have 
created standard profiles for each epoch of data, and fit- 
ted sets of 1,2, or 3 (as necessary) gaussians to each of 
the pulse and interpulse (seperately). All fitting was done 
using the gaussian fitting package bfit [3*]. After experi- 
menting with several different methods of aligning the 
gaussian templates, we decided that the smoothest align- 
ment was achieved by keeping the phase of the tallest 
gaussian constant. This also seems to produce a fairly 
monotonic behaviour in the phase of the interpulse. The 
full collection of modelled profiles is shown in Figure 1 . 
These were used as standard profiles for timing. 

TEMPLATE EVOLUTION 

The changing profile shape poses a problem in the cre- 
ation of times-of-arrival (TOAs) for the system, as the 
changing profile shape allows for ambiguity in the fidu- 
cial point of the profile, which translates into additional 
timing noise for the system. Figure 2 shows the changing 
fiducial point of the gaussian profiles over time. In order 
to reduce this source of timing noise, we have decided 
to 1) use a series of standard profiles developed from the 
well-modelled epochs of ASP, GASP and WAPP data, 
and 2) align the Gaussian templates as described in the 
previous section. 



ORBITAL ABERRATION 

We obtained GBT time over four days between October 4 
and 12, 2006, in an intensive campaign that has allowed 
us to study the effect of orbital abeiTation without the 
effects of the long-term profile evolution. Over such a 
small time span, the secular evolution of the pulse profile 
is negligible, and so we can be confident that any changes 
over orbital phase bins (as shown in Figure 3) are due 
solely to the change in observed emission cone direction 
as the pulsar moves at a relativistic orbital velocity. 

TIMING 

We present timing using L-band data from telescopes at 
Green Bank, Arecibo and Nancay. For each observation, 
we use the Gaussian template nearest in time. A total of 
1676 TOAs were fitted using the pulsar timing package 
TEMP02 Even with the multiple templates, there 
was still a large amount of timing noise. We removed this 
using the 'fitwaves' feature of TEMP02 Jsj], which fit out 
13 harmonically -related sinusoids as illustrated in Figure 
4a, resulting in the residuals presented in Figure 4b. Ta- 
ble 1 shows the preliminary orbital and mass parameters 
resulting from this fit. More work is needed on the tim- 
ing of this system since, for example, we still need to 
coiTect for frequency-dependent profile differences that 
are apparent in some of the Arecibo data. From our tim- 
ing analysis, we obtain an improved companion mass of 
1.365±0.018 Mq, allowing us to more confidently as- 
sume that the companion is a neutron star. 
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FIGURE 1. The fits of the pulse and interpulse of 

J1906+0746, based on ASP (Arecibo; coherent dedispersion), 
GASP (GBT; coherent dedispersion) and WAPP (Arecibo; 
spectrometer) observations, each modelled with several gaus- 
sian components. The pulses are spaced vertically based on 
observation date, and the interpulse is magnified by a factor 
of ten relative to the main pulse. Different machine dispersion 
methods account for some but not all of the apparent interpulse 
differences in nearby observations; the re maining differences 
are under investigation. The vertical line for the main pulse 
(left) illustrates the chosen alignment, whereas the vertical line 
for the interpulse demonstrates the marlced phase shift of the 
interpulse relative to the main pulse (see proceedings by Desvi- 
gnes et al.) 
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FIGURE 2. Phase offset of the profile at each of 29 epochs 
of data. Offsets computed by comparing one epoch of data with 
each of 29 gaussian-modelled standard profiles (as described in 
'Gaussian Fitting and Profile Evolution' section). A regression 
line is plotted through the points, and although there is a fair 
amount of scatter in the data, it roughly follows a Unear trend. 
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FIGURE 3. J 1906+0746 pulse and interpulse for 8 orbital 
phase bins for two days of data are shown in the two leftmost 
sets of plots - a third day will be incorporated once it has 
been recalibrated. Note that the interpulse has been magnified 
by a factor of 8 relative to the main pulse. The two columns 
of plots on the right show the difference between the profile 
in each orbital bin and an average profile summed over the 
two days of data. The difference profiles show no evidence of 
profile changes on an orbital timescale within the noise. We 
will use these limits along with the long-term profile evolution 
to constrain the emission beam 
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FIGURE 4. a) The timing noise required 13 harmonically- 
related 'fitwaves' sinusoids to be fit out using TEMP02, with 
a fundamental period of 2.322 years. The resulting variations 
attributed to timing noise are plotted in the top panel. In the bot- 
tom panel, we have subtracted from these the best-fit parabola 
to demonstrate the large short-term variations. Based on the 
amplitude of the parabola, we estimate that the use of fitwaves 
has affected the fit period derivative at the level of about 2.5%. 
b) Timing residuals using TEMP02 and a relativistic timing 
model (Damour & Deruelle 1986) | 6]. Residuals are black: 
Nancay; dark grey: Arecibo (WAPP); and light grey: GBT 
(GASP). 
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FIGURE 5. This mass-mass diagram displays constraints 
only from the advance of periastron to and the time dila- 
tion/gravitational redshift parameter y. Adding the orbital pe- 
riod derivative to the fit changes the fit only marginally, and the 
fit does not yet appear to be robust. The black dot represents 
the masses derived assuming general relativity. 

TABLE 1. Preliminary Mass and Orbital Parameters for 
J 1906+0746 

(h 7.58331(95) 

r 0.0004930(84) 

Ml 1.248(18) Mq 

M2 1.365(18) Mq 
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